gusucode.com > MATLAB与控制系统仿真实践 源代码 > MATLAB与控制系统仿真实践 源代码/第十四章/eg14_9.m

    num1=5;
den1=conv([1 2],[1 5 0]);
G1=tf(num1,den1);
margin(G1)
%%%%%%%
figure(2)
num2=200;
den2=conv([1 2],[1 5 0]);
G2=tf(num2,den2);
margin(G2)

%分别得到系统幅值稳定裕度Gm，系统相位稳定裕度Pm，
%得到系统相角穿越频率wcg，系统幅值穿越频率wcp
[Gm,Pm,wcg,wcp]=margin(G1)
GmdB=20*log10(Gm);                     %系统幅值稳定裕度Gm，可以通过20*log10(Gm)转换为GmdB
[Gm,Pm,wcg,wcp]


%%%%%%%%%%%%%%%
bode(G1)